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Motivation 


•  Feedback  cycle  between  liquid  rocket  engine  (LRE)  combustion  chamber  pressure 
perturbations  and  unsteady  combustion12 

•  Large  amplitude  fluctuations  in  pressure  and  combustion  heat  release  rates  =^> 
combustion  instability 


Flow 

^■1 

Heat  Release 

Perturbation 

Oscillation 

Pressure 

Perturbation 


'Harrje,  D.T.,  and  Reardon,  F.H..  Scientific  and  Technical  Information  Office ,  National  Aeronautics  and  Space 
Administration,  NASA  SP- 194,  1972. 

2Schadow,  K.C.,  Gutmark,  E.,  Parr,  T.P.,  Parr,  D.M.,  Wilson,  K.J.,  and  Crump,  J.H..  19th  AIAA  Fluid 
Dynamics,  Plasma  Dynamics  and  Lasers  Conference,  AIAA  1987-1326 
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Motivation  (cont’d) 
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Injector  head 
molten  after  a 
combustion 
instability 
event 


F>  1  injector  heads, 
before  and  after  a 
bout  with  combustion 
instability.  Notice  how 
the  surface  of  the 
burned  Injector  plate 
is  peeling  away. 

right) 
Before,  *f»  - 
After,  nght 


wmMw-M 


Courtesy:  U.S.  Rocket  and  Space  Center,  Huntsville,  AL 
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Objective 


Impose  external  acoustic  perturbations,  and  examine  the  response  and  stability 
characteristic  of  shear-coaxial  injector  flow  to  pressure  perturbation 


Acoustic/Pressure 

Perturbation 


Shear-Coaxial 
Injector  Flow 


•  Investigate  influence  of  injector  geometry  on  flow  response  to  external  pressure 
perturbation 

•  Vary  the  outer- to-inner  jet  momentum  flux  ratio,  J,  under  subcritical  chamber  pressure 
condition,  i.e.,  reduced  pressure  Pr  =  0.44 


=  J^hambeL_  =  493psi  (3.4MPa) 


Pr  = 


p 

critical,  N2 


Characterize  mixing  using  dark-core  length  measurements 

Apply  proper  orthogonal  decomposition  of  high-speed  image  pixel  intensity 
fluctuations  to  extract  spatial  and  temporal  characteristics  of  prevalent  coherent  flow 
structures 
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Schematic  of  Experimental  Facility 
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Image  of  Experimental  Facility 
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Additional  Geometric  Configurations 


LAR-thin 


Similar 


Similar  post 
thickness 


SAR-thin 


8^ 

Ml 

I.D. 

sp- 
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LAR-thick 


I.D. 
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SAR-thick 


All  dimensions  in  mm 


|  Injector 

D, 

t/Dl 

AJA, 

l/Dl 

LAR-thick 

0.51 
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2.42 

3.18 

1.05 

12.9 

0.5 

SAR-thin 

1.40 

1.65 

2.44 

3.94 

0.09 

1.6 

0.5 

SAR-thick 

1.47 

3.96 

4.70 

6.35 

0.84 

2.9 

-0.1 

|  LAR-thin 

0.70  0.89 

2.44 

3.94 

0.13 

10.6  -0.2 

SAR,  LAR  —>  Small,  Large  Area  Ratio 
Thick,  Thin  —>  Post  thickness 


Davis,  Rodriguez,  Leyva  et  al. 

Rodriguez,  Graham  et  al.  (///),  =  0) 

Present  study  (-0.11  Dt  recess) 
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] 

Present  study  (-0.21  Z>j  recess) 
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Acoustic  Field  Set-Up:  Pressure  Antinode 


•  Pressure  antinode  (PAN)  -  condition  of  maximum  pressure  perturbation  in  the 
acoustic  field 

•  Piezo-sirens  forced  in-phase 

•  Superposition  of  quasi- ID  acoustic  waves  traveling  in  opposite  directions  =^>  PAN  at 
the  jet  location  (geometric  center  of  test  section) 


Dark-Core  Length  Measurement 


•  First  raw  grayscale  images  were  converted  to  binary  images 

•  A  contour  was  drawn  around  the  “dark-column”  in  the  binary  image 

•  Axial  length  of  the  dark-column  measured  and  defined  as  the  Dark-Core  Length,  L 


Grayscale  (8-bit)  Binary  Image  Raw  Image  w/ 
Raw  Image  Contour 


Baseline  Dark-Core  Lengths 
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Proper  Orthogonal  Decomposition 


UU  I  PMCBVM  ■  Mvli^  i 


•  Proper  Orthogonal  Decomposition  (POD)  or  Principal  Component  Analysis 
(PCA)  was  used  for  extracting  dominant  dynamical  processes  embedded  in  high¬ 
speed  images. 

•  A  time-resolved  set  of  images  A(x,t )  can  be  represented  as  a  linear  combination 
of  orthonormal  basis  functions  (j)k  (aka  proper  orthogonal  modes)1’2’3  : 

M 

A{x,t)  =  YJak(t)<pk{x) 

k= 1 

where  ak(t)  are  time  dependent  orthonormal  amplitude  coefficients  and  M  is  the 
number  of  modes 

•  Main  idea:  POD  modal  amplitudes  capture  the  maximum  possible  “energy”  in 
an  average  sense4,  i.e., 

k  k 

where  bk(t)  are  the  temporal  coefficients  of  a  decomposition  with  respect  to  an 
arbitrary  orthonormal  basis  y/k. 

1  Chatterjee,  A.  Current  Science ,  Vol.  78,  No.  7  (2000) 

2  Arienti,  M,  and  Soteriou,  M.C..  Phys.  Fluids  21,  1 12104  (2009) 

3  Narayanan, V.,  Lightfoot,  M.D.A,  Schumaker,  S.A.,  Danczyk,  S.A.,  and  Eilers,  B..  ILASS  Americas,  2011 

4  Berkooz,  G.,  Holmes,  P.,  and  Lumley,  J.L..  Annu.  Rev.  Fluid Mech.  25.  539  (1993) 
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Construction  of  Data  Set 
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•  First,  form  a  row  vector  consisting  of  all  pixel  intensity  values  of  each 
snapshot  image  (with  resolution  of  n  rows  by  m  columns)  in  order  of  increasing 
columns,  then  increasing  rows 


Image  Frame 


columns 


I 

n 

rows 

i 


•  Then,  combine  all  such  row 
vectors  for  N  sequences  of  image 
frames  resulting  in  a  matrix  A 
consisting  of  N  rows  by  P  =  n  x  m 
columns  of  intensity  values. 


\7 


•4 


_ P  =  n  x  m 

pixel  intensities 

N  time 
steps 
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Orthogonal  Decomposition  Technique 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 

•  SVD  preferred  since 

1 .  Applicable  to  non-square  matrices  (most  likely  the  case) 

2.  Decomposition  matrices  are  orthogonal 

3.  Subroutine  readily  available  in  MATLAB® 


•  Subtracted  temporal  mean  of  A  matrix  of  intensity  fluctuations  A 


•  Applied  SVD 


^A  =  USVT  =qvt 

Columns  of  Q  ~  ak(t )  contain  temporal  information 
Columns  of  V  ~  (j)k{x)  contain  spatial  information 


NxN 

Orthogonal  Matrix  of  Left 
Singular  Column  Vectors  of  A 


N*P 

Orthogonal  Matrix  of  Right  Singular 
Column  Vectors  of  A  proper 

orthogonal  modes  (POM) 


NxN 

Diagonal  Matrix  of 
Singular  Values 
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Results  -  Subcritical  Baseline  at  Low  J 


I4i*6  i 


LAR,  Pr  =  0.44,  J=  0.5 

Snapshot  Average 
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Identified  with  Characteristic 
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Amplitude  information 
contained  in  singular  values 
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of  Temporal 
Coefficients  of 
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Results  -  Subcritical  PAN  at  Low  J 
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•  LAR,  Pr  =  0.44,  J  =  0.5,  forcing  Frequency,^  =  3.14  kHz 


Symmetric  Structures 

Identified  with  Characteristic 
Frequency  at  /F 


Cross-Power  Spectral  Density  (CPSD) 


UU  IBHOHH  ■  frtrtie*  i  Mu vfft 


CPSD  yields  the  FFT  of  the  cross-correlation  of  the  temporal  coefficients 

cov(ar,ai;)  .  7 

CPSD  =  X 


£=0 


cr  .  a 


a 


a , 


•  Magnitude  and  phase  plots  used  to  determine  existence  of  propagating  structures 

LAR,  Pr  =  0.44,  J  =  0.5  PAN 

Baseline  (/F  =  3. 14  kHz) 

Baseline  Spectrum 
Completely  Removed  in 
PAN  Forced  Spectrum  |  io‘ 


-90°  Phase  Difference 
Confirmed  Downstream 
Propagating  Flow 
Structures 
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Sample  Animation  -  PAN  (fF  =  3.14  kHz) 


•  LAR  Pr  =  0.44,  .7=0.5 


Raw  Images 
Time  =  0.00  ms 


POMs  1  and  2  POM  1  POM  2 

Time  =  0.00  ms  Time  =  0.00  ms  Time  =  0.00  ms 


t 

Superposition  of  POMs  1  and  2  Resulted  in  Downstream  Propagating  Structures 
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Results  -  LAR-thin 


Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities  for  all  J 
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Non-Dimensional 
Spatial  Growth  Rate 


Previous  Works  on  Jet  Instability 


Michalke  and  Hermann  (1982)  did  linear,  inviscid  instability  analysis  of  a  circular  jet 
with  coflow 

Showed  that  with  increasing  coflow  velocity,  UrX) 

—  Helical  disturbances  more  unstable  than  axisymmetric  ones  farther  downstream  of  exit 

—  Jet  flow  becomes  less  unstable,  but  spectrum  of  spatial  growth  rate  becomes  broader  and  the  peak  shifts 
to  higher  frequencies 


Increasing  Downstream  Distance 
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O-l  0-2  0-3  0-4  0-5  0-6 

(W  _ _ 

sv 


Helical 

Axisymmetric 


->  Non-Dimensional  Frequency 


•  Dahm  et  al.  (1992),  Wicker  and  Eaton  (1994)  conducted  experimental  investigation  of 
large-scale  vortex  structures  in  the  near  field  of  coaxial  jets 


Dahm  et  al ,  JFM  1 992 


For  outer-to-inner  jet  velocity  ratios  greater  than  one,  found  that  coherent  structures  in  the 
layer  dominate  those  in  the  inner  shear  layer 

At  large  axial  distances,  shear-layer  vortices  exhibit  helical  structures 
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outer  shear 


CPSD  Magnitude 


Results  -  SAR-thin 


Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 


J=  2.0 


/=  17 


Baseline 


PAN 


Baseline 


PAN 


(/f  =  2.97  kHz)  (/p  =  2.90  kHz) 


Results  -  SAR-thick 


Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 
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Summary  and  Conclusions 


•  Examined  the  effect  of  different  exit  geometries  on  the  mixing  characteristics  as  well 
as  the  behavior  of  flow  disturbance  structures  with  and  without  transverse  acoustic 
forcing  —  mixing  not  only  depends  on  momentum  flux  ratio ! 

•  Proper  orthogonal  decomposition  of  high-speed  image  intensity  fluctuation  data 
revealed  key  spatial  and  temporal  characteristics  of  flow  structures 

•  Low  /,  SAR  injector  flows,  had  significantly  lower  LB/Dl  than  the  large  area  ratio 
flows 

•  Low  JLAR-thin  injector  flow  showed  a  strong  response  at  the  PAN  forcing  fre¬ 
quency  while  the  high  J  appeared  less  responsive  and  retained  the  baseline  flow 
spectral  characteristic 

•  SAR- thin  injector  flow  showed  strong  response  at  the  PAN  forcing  regardless  of  J 

•  Low  J  SAR- thick  injector  flow  showed  no  response  at  the  forcing  frequency,  while  the 
high  J  flow  did 

•  Operated  at  high  enough  J,  LAR  injector  flows  less  vulnerable  to  external  pressure 
disturbances 


22 
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Data  Summary  Table 


Injector 

J 

R 

Tch 

K 

Pch 
kg/m 3 

Pch 

MPa 

T 

1  0 

K 

ni0 

mg/s 

Po 

kg/m 

ua 

m/s 

Re0 

xlO4 

K 

rh t 
mg/s 

Pt  3 
kg/m 

u, 

m/s 

Re, 

xlO4 

SAR-tliin 

2.0 

6.9 

246 

21 

1.49 

195 

450 

27 

6.6 

1.1 

109 
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630 

0.96 

1.5 

SAR-tliin 

17 

20 

217 

24 

1.48 

194 

1300 

27 

19 

3.1 

108 

925 

638 

0.95 

1.4 

LAR-thin 

2.1 

7.4 

220 
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205 
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2.1 

7.3 
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25 
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24 

15 
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2130 
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3.3 
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21 

23 
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27 
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4385 

25 

35 
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1658 
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25 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


